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Abstract 

In this paper, we review the production cross section for charged and neutral Higgs 
bosons pairs in e+e~ collisions beyond the tree level, in the framework of the 
Minimal Supersymmetric Standard Model (MSSM). A complete list of formulas for all 
electroweak contributions at the one loop level is given. A numerical code has been 
developed in order to compute them accurately and, in turn, to compare the MSSM 
Higgs bosons pair production cross sections at tree level and at the one loop level. 
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1 Introduction 



One major task of a future e~^e~ linear collider will be the exploration of the Higgs sector, 
in the Standard Model and beyond it. Despite its success, the Standard Model suffers 
from the appearance of quadratically divergent contributions to the Higgs boson mass. 
However, this problem is solved by Supersymmetry. The theoretical framework of our study 
is the minimal supersymmetric extension of the Standard Model (MSSM), where the Higgs 
spectrum consists of three unphysical Goldstone modes (G"*", G~ and G^) as well as five 
physical states. Two of these are charged {H~^ and H^) and, among the three neutral Higgs 
bosons, two are CP-even states, and H^, and one is CP-odd, A^. Several other aspects 
of the MSSM Higgs boson phenomenology are reviewed in [1]. 

The processes e~^e~ H~^H~, H^A'^, h^A^ that will be observable at future e'^e~ linear 
colliders, such as ILC and/or CLIC, are among the best places where one can accurately 
check the Higgs structure, see references [2] to [8] for details. At tree level, e+ and e~ 
annihilate through a photon and a Z boson in the case of H^H^ production, and through 
only a Z boson in the case of H^A^ and hPA^ production. At this level, the amplitudes 
depend on the masses of the Higgs bosons and on the mixing angle a. At the one loop 
level, most of the MSSM parameter space is involved through self- energies, triangle and 
box diagrams. In a previous paper [9] it was shown that, at high energy, at the leading and 
sub-leading (Sudakov) logarithmic orders, a great simplification occurs. The gauge and the 
SUSY structures of these processes reflect directly in the coefficients of the quadratic and 
linear logarithmic terms. In this high energy range, they depend only on a few parameters 
(the Standard Model inputs, the angles a and /9, as well as the SUSY scale Msusy)- 
The next step is to study more deeply the SUSY structure by looking at sub-sub-leading 
effects. First, one should determine the energy range in which the above Sudakov limit 
is an acceptable approximation and can be accurately tested. Then, one can study the 
effects of the successive sub-sub- leading terms (constants, m^/s, etc) and classify the 
various parameters which control each of them. We should then estimate the accuracy 
at which these parameters can be measured. For these purposes, we have developed a 
code allowing to compute numerically the complete electroweak one loop contributions 
to the pair production cross section of MSSM charged and neutral Higgs bosons in e+e~ 
collisions. The purpose of the present paper is to write in an exphcit fashion all details of 
the electroweak one loop contributions that are computed by this numerical code. 

In Section 2, we review the tree level MSSM Higgs sector and we calculate the production 
cross section for H'^H" , H^A^ and hPA^ pairs in e^e" coUisions at tree level. In the rest of 
the paper, we focus on the various one loop terms. The contributions of the initial vertices 
and of self-energy are given in Section 3, the intermediate gauge boson self-energies 
are discussed in Section 4, the contributions of final vertices and of Higgs self-energies are 
calculated in Section 5, and the effect of box diagrams are presented in Section 6. Finally, 
a summary and some outlooks (in particular a more detailed description of our numerical 
code) are given in Section 7. 
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2 Tree level calculations 



2.1 Tree level structure of the MSSM Higgs sector 

In the MSSM, two complex scalar Higgs doublets are responsible for the breaking of the 
electroweak symmetry: 

They have opposite hypercharge (Yi = — 1 and Y2 = +1) and their vacuum expectation 
values are respectively vi and V2. After diagonalization, one obtains the following states: 




cos a sm a 

— sin a cos a 

cos P sin (3 

— sin /? cos /? 

cos /? sin (3 

— sin (3 cos (3 




(2) 

(3) 
(4) 



Here, G"*", G~ and describe three unphysical Goldstone modes. The five physical states 
are two charged bosons {H^ and H~), two neutral scalar bosons with CP = +1 and 
and one pseudoscalar neutral boson with CP = —1 {A^). 

The quadratic part of the Higgs potential, which contains the soft breaking masses and the 
gauge couplings, depends on two independent parameters, which are usually chosen as the 
mass Ma of the boson and the ratio between the vacuum expectation values tan /3 = V2/V1. 
The masses of the other physical states are expressed as follows: 

M| = Ml + M^, (5) 
M'HO,hO = ^ (m1 + M| ± ^(Ml + M|)2 - 4M1M| cos^ 2/3^ . (6) 
As for the mixing angle between and it is given by: 

tan2Q; = tan 2/3 x ^ tlf , - - < a < 0. (7) 

Note that these results are only valid at tree level and they become slightly different when 
one includes radiative corrections. 



2.2 Production cross section at tree level 

In e+e~ collisions, charged Higgs bosons are pair produced through virtual photon and 
Z boson exchange (and in top decays if Mh is small enough). As for the neutral Higgs 
bosons, they can be produced through several mechanisms: 1^1^ and ZZ fusion processes, 
Higgsstrahlung or pair production. In this paper, we only focus on this latter process (note 
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that CP conservation forbids virtual photon exchange) . The Feynman diagrams of interest 
are shown in Figure 1. More details about the various production mechanisms and decay 
modes of MSSM Higgs bosons can be found in [10]. Here, we only focus on the total pair 
production cross sections and we ignore the different contributions of the decay channels. 









e" 


* 










e 







Figure 1: Feynman diagrams for the pair production of MSSM charged and neutral Higgs 
bosons in e+e" collisions. 

The tree level production cross section can be easily derived using the Feynman rules. If 
sw = sin^vK) c-w = cos6';y and r] = q^/{q^ — M^), then the Born amplitudes af *^" are: 

• for charged Higgs bosons: 

afyiH^H-) = 1 - ^-^r=^V9L,n (8) 

• for neutral Higgs bosons: 

afy(H'Ayh'A') = -^L_rjg,,n [Z„,] (9) 

where = 2s^ - 1, = and [Zab] = [-sin(/3 - a); cos(/3 - a)] for H°A° and h°A° 
final states, respectively. 



3 



In this paper, we use the following renormalization for the amplitude: 



v{e+){i>){aLPL + aRPR)u{e-), Pl,r = (10) 
q z 



The differential tree level cross sections are then given by: 

,j Born 2 a3 

aa^ = X (1 - cos^ e) X \all^\\ (11) 

rfcos^ 8g2 ^ I \ L,R \ \ ) 

Here, j3H is the velocity of the outgoing Higgs bosons. If Afi and M2 are the masses of the 
two outgoing Higgs bosons, then (3h{Mi, M2) is defined by: 

2\p\ 1 



yS s 



\ ^ M2 + Mi\ / M2 + MA ( M^-MA / M^-MA 



After integration over cos 6', one gets: 

• for charged Higgs bosons: 

(, 4Miy/\ / 2cW , c-/(4 + 4) \ .... 

''^^^-^48^i^"^J H ^^^M|77 + (1 - M|/.)2 J 

_ -l+4s2^ -1 , -l + 2s^ 

with Cv = — , CA = , Cy = — . 

ASwCw ASwCw ZSwCw 

• for neutral Higgs bosons: 



Born _ ^ ^ ^"W ^"W ' ^ \ ^ \7 \^ ^ r^jf \--^n-/n", ^-^ji) 

'^HOAO/HOAO - ^ X 32,4^,4^ ) >< >< (1 _ M|/.)2 " ^^^^ 

In the decouphng limit (M4 ^ M^- and Ma — M^"), cos(/5 — a) — * and e+e~ ^ /I'^A'^ is 
strongly suppressed, i.e. only the H^A^ pairs can be produced in e~^e~ collisions, with a tree 
level cross section given by: 

Figure 2 shows the pair production cross section for the MSSM charged and neutral Higgs 
bosons in e'^e~ collisions, at tree level, as a function of Ma and for various values of the 
centre-of-mass energy. 
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Figure 2: Tree level pair production cross section for charged (top) and neutral (bottom) 
Higgs bosons in e+e" collisions, as a function of the MSSM parameter Ma and for various 
centre-of-mass energies y/s. For simplicity, we performed our calculations in the decoupling 
limit, where H, and are almost degenerate in mass. 
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2.3 Complete amplitude for calculations at the one loop level 

The analytical expressions of all electrowcak one loop contributions to the MSSM Higgs 
bosons pair production cross sections are given in the following: e"^ self-energy and initial 
vertices in Section 3, 7 and Z self-energies with counter terms in Section 4, final vertices and 
Higgs self-energy in Section 5, and box diagrams in Section 6. The complete renormalized 
amplitude used to calculate the cross section is the sum of Born and one loop terms: 

-> Higgs pair) 
Higgs pair) 

■ Higgs pair) -|- A'^^{e'^e~ Higgs pair) 
Higgs pair) 

Higgs pair). (16) 

The notations used in our calculations, in particular when writing vertices in terms of real 
coupling constants, are described in Appendix A. In the following, we use a formalism which 
involves Passarino-Veltman functions, see Appendix B for details. 



A(e+e- ^ Higgs pair) = A^°™(e+e- 

+ A'"'%e+e- - 



3 Contribution of initial vertices and self-energy 

The amplitudes corresponding to initial triangles and self-energy are: 

• for charged Higgs bosons: 

alniH^H-) = i— + x (17) 

e iswCw e 

• for neutral Higgs bosons: 

where we write r*"'^ = r^^'^P^ + ^r'^Pr for ^ = 7 or Z. 

pin, 7 g^j^^ pin,z ^jjg same in the charged and neutral sectors, since they only depend on 
the initial state. They are obtained by summing various contributions: 

pm,y ^ r^+^_{WiyW)-r^+^-ipinch) 

+ r^+e- (lyWu) + r^+e- (eZe) + F^^,- (e^e) 

+ r,V(™i^L) + rr+e-(ex°e) 

+ r,V(e.s.e). (19) 

The particles inside the initial triangle have internal masses mi, 1712 and 1713. They are 
ordered clockwise, mi being the mass of the particle just after the junction involving the 
momentum q. 
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1) The contribution of the WuW triangle is: 

where Cww = -l2C2i{WuW) + 2-2q^ [CoiWvW) + CuiWvW) + C23(M/z/H^)]. 

2) In e"'"e~ annihilations, VFVT contributions arise in the photon and Z self-energies, as well 
as in the triangles connecting the photon and the Z boson to the initial e^e^ pair or to 
the final Higgs pair. Therefore, it is convenient to extract a certain part (so-called pinch) 
from such a triangle with two W lines (in our case the WvW triangle) and then to put 
it inside the photon and Z self-energies contributions, in order to have universal charge 
renormalization [11]: 

.r^<^^,) - - - { \!:z forV = z ■ (21) 

3) As for the vWv triangle, since neutrinos do not couple to photons, one has: 

Tl,^.{vWv)^Q (22) 

while, for the Z boson, one obtains: 

where Cw = AC2a{vWv) - 2 + 2?^ ]^Cii{vWv) + C2^{uWv)]. 

4) The contribution of the eZe triangle is: 

n.e- ieZe) = T^^Cz [alPL + glPn] (24) 

or 

rf,,_ (eZe) = - :^^^^C, [glP, + glPn] (25) 
where Cz = ^C2i{eZe) -2 + 2q^ [Cn{eZe) + C23{eZe)]. 

5) The contribution of the 676 triangle is: 

r:.e-(e7e) = ^C,[P^ + P«] (26) 

or 

rf,,_ (e7e) = - -^^C, [g^PL + gRPn] (27) 

QTTSwCw 

2 



where = 4C24(e7e) - 2 + 2g^ [Cii(e7e) + C23{e-fe)]. 
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6) The contribution of the Xj^LXi triangles is: 



e+e 



iXii>LXi) 



E 



or 



^e+e-{Xii>LXi) 



where = C2A{Xj^LXi) ~ 4 + y [^^'i^Xj^LXi) + C2^{Xj^LXi)\ and Cl^ = Co{Xjj>LXi, 
7) As for the x°ex° triangles, since neutralinos do not couple to photons, one has: 

r:+e-(x;ex°) = o 



while, for the Z boson, one obtains: 



^3 



X 



K'iPr. + K'iPR 



by defining K^^ and K^^ as follows: 



(Z^/sw + Z^*cw){Z^iSw + Zgcw) 



X 



2(^i;^3? - ^£^4r)C'I + M,oM,o(Z3^*Z3^ - ZZ^Z^^Ci^' 



where = C^.ix'.exi) + J [Ci2ix',ex^) + C^six^xf)] and = Co(x?e~x°). 
8) As for the i'LXi^L triangles, since sneutrinos do not couple to photons, one has: 

rj+e-(z>LXi'>L) = 

while, for the Z boson, one obtains: 



{^LXi^L) 



E 



where C24 = C2A{i^LXi^L)- 
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9) The contribution of the e triangles is: 



or 



where (7^4 = C24(ex°e). 



10) The electron self-energy (e.s.e) contributions are obtained as follows: 
or 

rf+,_ (e.s.e) = [SlOlPl + W^'fl] 

where the following loops are taken into account: (Wv), (Ze), {'ye), (x^^), (x°e). 

For the (Wu) loop, one has: 



5r(W^i/) = 0. 



For the {Ze) loop, one has: 



*.(^e)^-^(A(2e.0) + l). 



For the (7e) loop, one has: 



?>L{ie) = 5«(7e) = - ^ (Bi(7e, 0) + . 



For each (Xj^^) loop, one has: 



For each (Xi^) loop, one has: 



= -^^\Z^^sw + Z^,cw\'B,{xreL,0), 
Mx°e~) = --^|Z,tri?,(x°e~^,0). 
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4 Contribution of 7 and Z self-energies 
4.1 Definition of gauge self-energy functions 

The on-shcll renormalization procedure [12, 13] that aUows fuU determination of the MSSM 
Higgs sector at one loop, as well as of the corresponding counter terms, makes use of several 
gauge self-energy functions, which are detailed in the following of this section [13, 14, 15]. 
Let us first define several useful expressions: 

A//"2 _|_ A//"2 ^2 

PV\{XY,q^) = ^ ; ^ - 2i?22(Xy, q^)-%r- \B,{XY, q') + B,,{XY, q')] , (49) 

2 L -I 

PV2{XY,q^) = 10S22(Xy,g2) + (4g2 + M2 +M2)So(Xy,g2) 



2' 



+ A{Ml) + A{Ml) - 2 I^Ml- + - |- j , (50) 

PV3(xy.,^) = 2BM^)-^^m±Mm^(^^^zmB„(xy,f). (51) 

a) Photon self-energies: 

The photon self-energy is defined as: 

^77(5^) = ^77(5^) + A^j{pinch). (52) 

The pinch term is given by: 

A^^ipinch) = -^q^B^iWW, q^). (53) 

TT 

The self-energy term without pinch S^^(g^) is the sum of various contributions: 

i:,,{q^) = E^^(g+H) + E^^(//) 

+ E^^(xx) + S^7(/7)- (54) 

The contribution of the gauge and Higgs sectors is: 

E^^(g+H) {2B22{HH,q^)-A{Ml)+&B22{WW,q^) 

-3^(M^) + 2q^Bo{WW, ?') + j}- (55) 
The contribution of the fermion pairs is: 

S77(//) = E^^^^{^^l(//,?^) + M|So(//,g^)}. (56) 
/ ^ 

The contribution of the chargino pairs is: 

S77(XX) = T.'^{Py^ix^Xi.Q') + MlMxiX^.Q% (57) 

i 

10 



The contribution of the sfermion pairs is: 



(5J 



j=l,2 



Here, fi, /2 account for fi, in the case of unmixed sfermions, or for the physical states 
obtained after mixing (i.e. ti, ^2, ^2 in the case of third generation squarks). The 
couphng between a photon and a sfermion pair is the same with and without mixing. 

b) Z self-energies: 



The Z boson self-energy is defined as: 

Azz{q^) = T.zz{q^)+ Azz{pinch). 

The pinch term is given by: 

Azz{pinch) = _^£!^(52 _ M|)5o(W, q'). 



2^ 



(59) 



(60) 



The self-energy term without pinch Tizz{(f) is the sum of various contributions: 

Szz(g') = S^z(g+H) + E^^(//) 

+ Szz(xx) + Ezz(xV) + Szz(/7)- 

The contribution of the gauge and Higgs sectors is: 



(61) 



-'ZZ 



;g+H) 



{\ [A{Ml.) + A{Mlo) + A{Ml) + A{M: 



+ sin2(/3 - a) [MlBoiZh', q') - B^^{Zh\ q') - B22{A'H', q') 
+ cos\P - a) iMlBoiZH', q') - B^,{ZH\ q') - B^,{A%' , q') 



--cos\2ew) [2B^^{HH,q^)-A{Ml) 



8c^ + cos'(2^w^)] B22{WW, q^) 

AcW + cos{2ew)] BoiWW, q^) 



+- 



124 - 4cl. + 1 



w 



(62) 



The contribution of the fermion pairs is: 



Szz(//) = E 



{{g'yf + glf)PVl{Ml q') + {gl^ - g%)MfBo{ff, q')} (63) 



where gvf = T|(l - ^IQjls^,) and gAf = Tf. 
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The contribution of the chargino pairs is: 

2 



ij ij ' ij ij 



■\Z R* /nZ R 



The contribution of the neutralino pairs is: 



ij ij ' ij ij 



A.^ A. J 



The contribution of sfermion pairs is: 
• for unmixed sfermions: 



^ztiff) = - E ( ^ I {2i^22(//, q') - A{M})}, 



fL,R 



• with sfermion mixing (third generation squarks): 



27rs^c^ 



+2{T}^cj-si^QfyB,2{fifi,q') 

- 'c){Tl - sl^Qff + s)Qy^] A(M? ) 
+2(T^4-s2^Q^)2S22(/2/2,g') 

- - 'wQf)' + c}Q}sty] A{Ml)}. 



c) Mixed self-energies: 



The mixed self-energy is defined as: 

Ayz{q^) = ^^z{q^) + A^z (pinch). 

The pinch term is given by: 

A,z{pinch) = -^^{q^ - Ml)Bo{WW,q'). 

The self-energy term without pinch E^z(g^) is the sum of various contributions: 

E^z(?2) ^ E^z(g+H) + E^z(//) 
+ E^z(xx) + S^z(//)- 
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The contribution of the gauge and Higgs sectors can be expressed in several ways. Here, we 
choose the definition given in [15]: 



S,z(g+H) 



ftp 



47r 



SwCw 



B22{HH, q'') + B22iWW,q'') 



^2 _ „2 

SwCw 
Cw 



A{Ml)-rA{Mlr) 



+ 



Sw 



6A(M^) - 4M^ 



+2—B22{WW, q^) - 2swCwMlBo{WW, q^) 

Sw 

^^PV2{WW,q^)]. 
Sw ' 



The contribution of the fermion pairs is: 

aemNlQjgvf 



E 



2'KSwCw 



{PVl{ff,q') + M]Bo{ff,q')}. 



The contribution of the chargino pairs is: 

2 



ZR 
ij 



The contribution of the sfermion pairs is: 
• for unmixed sfermions: 

^0 



•sr\light 



iff) 



27r 



E ^iQi 



'Zff 



X 



fL,R 



• with sfermion mixing (third generation squarks): 



27lS\yCw 



E Qf {(The} - sl^Qf) [2B22{fJi, q') - A{Ml)] 



f=t,b 



(71) 



(72) 



(73) 



{2B22{ff,q')-A{M})}, (74) 



+{Tls} - s'^Qf) [2B22{f2f2, q') - A{Ml)\ }. (75) 



d) W self-energies: 



The self-energy term without pinch Yiww{(f) is the sum of various contributions: 

^ww{e,+^) + ^wwU'f') 
+ Eww{xx') + ^ww{ff'). 



(76) 
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The contribution of the gauge and Higgs sectors can be expressed in several ways. Here, we 
choose the definition given in [15]: 



{ - sin2(/3 - a) [B22(HH^, g") + ^22(1^/1°, g') 

- cos\p - a) [B22{Hh\ q^) + B22{WH\ q^) 
-B22{WZ,q^)- B22{HA\q^) 
+2sl,B22ilW, q^) + 2cl,B22{WZ, q") 

A{Ml,) + A{Ml,) + A{Ml) + A{Ml 

A(M^) + A{Ml 

sin2(/3 - a)5o(/i°W^, q^) + cos^{p - a)Bo{H^W, q^) 



+ 



+ 



Cw 

+ [3A(M2,) - + clr [3A(M|) - 2M 

-cl,PV2{ZW, q') - slrPV2{^W, q')}. 

The contribution of the fermion pairs is: 



1^ 47rs^ 



(77) 



(78) 



The contribution of the gaugino pairs is: 



'WW 



Up 



E { {or(^r* + oror*) py^mi g^) 
+ (^r^r* + (^r*^D %.^x?^o(x.x;, ?')}• (79) 



The contribution of the sfermion pairs is: 
• for unmixed sfermions: 



^TTS^^ (//') 



E^/ 



B22{ff,q') 



A{Mj)+A{Mj,] 



(80) 



with sfermion mixing (third generation squarks): 



^heavy 
-'WW 



iff) 



Sap 



{cfclB22{iibu q^) + c|s?522(ti62, g') 

+slclB22(t2h, ?') + slslB22(t2b2, ?') 



cfA(M^ 



c?A(M|J + 4A(Mfj]}. (81) 
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4.2 Charged Higgs sector 

For e^e~ H~^H~, the on-shell renormalization procedure leads to the following RG terms: 



R 



4s2 ^2 



^(1 - - (Jl.h) 



(82) 



with the corresponding counter terms: 



+ 
+ 



n^^(o) + 



cwMl 



X 



1 + vgL,R 



2 



X 



V 



4s2 ^2 



Ml 



X 



^PL + gnPR 



(83) 



Here, n^^(g'^) = — ^ (no pinch term) and n^^(O) is thus simply obtained as follows: 



n,,(o) 



77 



q2=0 



4.3 Neutral Higgs sector 



For e'^e 
terms: 



H^A^/h'^A^, the on-shell renormahzation procedure leads to the following RG 



al%{H^A^/h'A^) = - X 



rO AO luO AO 



T9l,r 



Azz{q') 



-A,z{q') 



(85) 



As for the counter terms, we only consider those corresponding to electroweak couplings and 
gauge boson masses here (the counter terms corresponding to H^A^ or hPA^ final states will 
be calculated in Section 5.3.9): 



- iv [Z, 



abl 



1 2Spj/ + 2Sy^r 



Pl + 



2r2 



R 



■'WW 



^zz{Ml) 



abl 



9l,r 

4s2 ^2 



n^^(o) + 



Ezz(M|) , 2 sw^^z{0) 
cwMl 



+ 



(86) 



5 Contribution of final vertices and Higgs self-energies 

5.1 Diagram structures for final triangles 

Several useful expressions are needed when estimating the contributions of the final vertices. 
The particles inside the final triangle have internal masses mi, 1712 and m^. They are 
ordered clockwise, mi being the mass of the particle just after the junction involving the 
momentum q. 
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Let Pfi and Pf2 (respectively Mi and M2) be the momenta (respectively the masses) 
two outgoing Higgs bosons (i.e. H^H~ or H^A^ or hPA^), then one has: 

Pf2 = Ml 

- (Mf + M|) 

^/l^/2 - ^ ■ 



a) Tril-type triangles: 



Cl — - + 6(Cooi — (^002) + PfiCiii — ^^2^222 

+ (2P/1P/2 — Pfi)Cu2 + {Pf2 — 2P/iP/2)Ci22 

+ 2 P/1P/2C21 — Pf2C22 + (-P/2 ~ P/1-F/2)C'23 ~ ^"24 

- {2Pf,Pf2 + P],){Cn-C,2). 



b) Tri2-type triangles: 



C2 = {SPf,Pf2 + 6Pl + 2Pf^)Co + {8PfiPf2 + 7Pl + P]2)Cn 

+ (P|i - P|2)Ci2 + (2P/1P/2 + 2P|JC21 + (2P/1P/2 + 2P|2)C22 

+ (4P/iP/2 + 2P|i + 2P|2)C23 + I2C24 - 2. 



c) Tri3-type triangles: 



— - + 6(Cooi — C'002) + PflClll — P/2C'222 + (2P/1P/2 — Pfi)Cii2 
+ {P% - 2PflPf2)Ci22 + P/1C2I - (2P/1P/2 + P/2)C22 

— 2PjiC23 — (fCi2 — 2C24 + -, 



C-i — Cii — C12, 



d) Tri4-type triangles: 



e) Tri5-type triangles: 



f) Tri6-type triangles: 



Cq — Co + Cii — Ci2- 



C4 — Ci2 — Cii — 2Co- 



C5 — Cii — C12 — Cq. 



Ce — Cii — C12. 
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5.2 Charged Higgs sector 

The amplitudes a{*'^ corresponding to final vertices with H^H~ states are: 



(9J 



2e 2swCw 2e 

Here, Tf''''^{H+H-) and T^''^^^ {H+R-) are obtained by summing the contributions of 
various triangles and of charged Higgs self-energy terms, as detailed in the following. 

For the photon, one has: 

+ ^^H+H- (3/) + ^^H+H- KXX^X) i- 

+ rW(4-leg) + r^+^_(i/.s.e). (99) 
For the Z boson, one has: 



r/->7(^+^-) = r^^^4ich) + r],^^_ (2) -r^+^_ (2, p^nc/.) 

tVh- m'x) + r^,^- (6ch) + r^,^_ (6/) 





^ H+H- 


(ich) + r|+^_(in) + r|+^_ 


(2)-r|+^_(2,pmc/i) 


+ 


^ H+H- 


(3/) + r|.^- (xx°x) + r|+^- (x°xx°) + r|,^- (4) 


+ 




(6ch) + r|+^_(6n) + r|+^_ 


(6/) 


+ 


i H+H- 


(4-leg) + r|+^-(i/.s.e). 


(100) 



5.2.1 Tril-type triangles 

The Tril-type triangles contribute to T^'^'^''^{H^H~) with: 



r^+^-(ich) 



87r2 



Ci(i/7i/)+ fV^l Ci{HZH) 



The Tril-type triangles contribute to rf''''^{H+H-) with: 



(Ich) 



( 



1 - 2,s-f 



Stt^ \ 2sm^cm/ 



A -24- 

y 2s\yC\y 



Ci{HZH) 



(101) 



(102) 



r|+^_(ln) =^^-^-^x {sin2(/?-a)[Ci(i/°V^^°)+Ci(^°W^i/' 

+ cos2(/3 - a) [Ci(/i°iyA°) +Ci(A°W^/i°)] } 

5.2.2 Tri2-type triangles 

With f ^ = i ^ ^ „ , „ , the contribution of the Tri2-type triangles is: 
1 cw/sw for y = Z ' b 



(2) 



X [sm^{l3-a)C2{WH^W) 

+ cos^{(3-a)C2{WhPW) 
+C2{WA^W)]. 



647r2s2^ 



(103) 



(104) 
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However, one must also take into account the pinch term: 



TH+H-(2,pinch) = 



X Bo{WW,q'). 



(105) 



w 



5.2.3 Tri3-type triangles 

The Tri3-type triangles contribute to both r-^^"'^(//+//-) and V^''''^{H+H~) with: 

1 



H+H- 



(3/) 



87r2 2sl,M^ 



E 



X 



{ [gvRfMjcoep + gvLfM} tan^ /?] C3(//7) 

- [^/viJ/'Mj, tan^/? + (7yi;,Mjcot2/5] C3(/7/') 

^-2M}M}{gvRf + C^(//7) 

-2M}M}{gvRr + ^vL/O C^(/7/') 

+MJ [^yi;M|cot2/3 + ^yi^/M|, tan^ /?] C^'(//7) 

-Mj, [^y^/.Mj, tan^ (5 + ^..^^^Mjcot^/?] C'i{f f f)} 



(106) 



-(xx°x) 



1 

87r2 



ijk 



CsiXkXiXj) 



(nVL R L* , nVR^L R* 
^kj '-Hji'-Hki ^ ^kj '-Hji'-Hki 



^R* 



^VRR L* 



C'siXkX^iXj) 



C'siXkXiXj) 
C':;iXkXUj)}, (107) 

where V — j or Z, and where (/, /') is either (g„, qa) or (i/^, £) for each fermion generation. 



^qVL L L* , (qVR R R* 
^kj '^Hji^Hki^ ^kj '-^Hji'^Hki 



In addition, Tri3-type triangles contribute to rf^"^'^{H'^H ) with the following term: 



rf,.H-{x'xx'' 



/-/^OL L L* I /nOR„R „R* 
^jk '-Hij'-Hik "I" ^jk '-Hij'-Hik 



Csi&xl) 



'^jk ^Hij^Hik ^ ^jk ^Hij^Hik 

n^L R L* . /oOi? L R* " 
^jk ^Hij^Hik "T *-^jfc '^-Hij^Hik 



C'six'jXiX'k) 

c',{x%xl) 

+M^oM^o [0'^^4,^4*, + Of C^'(x;x,X°,)}. (108) 



5.2.4 Tri4-type triangles 

The Tri4-type triangles contribute to V^'^^'^ [H^ H~) with the following term: 



H+H- 



(4) = 



1 (e^Mw{l 



167r2 



0„2 3 



X 



{gmHH cosiP - a) [CiiH^HZ) + C^iZHH^ 
+ghOHHsm{P-a) [C4{h''HZ) + C4{ZHh'')] }. 



(109) 
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5.2.5 Tri5-type triangles 

There is no contribution from Tri5-type triangles in the production of H^H~ pairs. 



5.2.6 Tri6-type triangles 

The Tri6-type triangles contribute to both r-^^"'T(i?+i?-) and V^''^'^{H+H-) with: 

.(6ch) [gyGG9\oGHC,{GA^G) 



^ H+H- 



+9VGG9H0GH^6iGH°G) + 9VGG9hOGH^&iGh^G) 

+9vhh91ohhC6{HH'H) + gyHH9loHHC,{Hh'^H)} (110) 



and 



The third generation squark contribution, with sfermion mixing, is: 

r5^r(6/) =i E {9vi.l,9m,l,9m,lM^ihh) 



;iii) 



ijk=l,2 



-9vutk9Ht~b, 9Hikb, ^6 (tibjtk) } . 



(112) 



The coupling of L-sfermions to the charged Higgs boson does not vanish like the fermion 
mass, so they also contribute to r^(6/). With (/, /') = {u, d), (c, s) or 3 x {ui, i), one has: 

r5S-*(6/) = ^ E - Svhh CM'Jl)}. (113) 



(//') 



In addition, Tri6-type triangles contribute to T^'^'^'^ {H'^H ) with the following term: 

87r2 \Asl^^cw) 



r|+^-(6n) 



{gnOGHsmiP - a) [CeiH'^GA^) + Ce{A'>GH' 
-9hOGH cos{(3 - a) [C6(/i°GA°) + Ce{AW)] }. (114) 

5.2.7 4- leg diagrams 

The 4-leg diagrams contribute to r-^*"''''(if+if~) with the following term: 



r^+^_(4-leg) 



87r2 



{2e' [Bo{H^,M]j)-B,{H^,MI] 



sin^(/9 - a) 
cos^{P - a) 



Ao2 



Bo{HZ,M]j)-Bi{HZ,Ml 
Bo{H^W, Ml) - B^{H^W, Ml) 
'Bo{h!^W, Ml) - Bi(h!^W, Ml) 
Bo{A'W, Ml) - B,{A^W, Ml)] }. 
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(115) 



The 4-leg diagrams contribute to T^'^'^'^[H^H ) with the following term: 



r|+^-(4-leg) 



w) 



Stt^ 2swcw 

sin^(/5 - a) 

'iSwCw 

cos^(/5 - a) 



Bo{H^,M],)-B,{H^,Mj,] 

'b^{hz,mI)-b^{hz,mI) 

Bq{H^W, Ml) - B^{H^W, Ml) 
Boih^W, Ml) - Bi{h^W, Ml) 



AswCw 

'BoiA'W, Ml) - B,{A^W, Ml)] } 



AswCw 

5.2.8 Charged Higgs self-energies 

The charged Higgs self-energies contribute to T^^^''^{H^H~) and T^^^''^ [H^ H~) with: 



(116) 



Vl+jj_{H.s.e) = 2e X 



(117) 



and 



Vl^H-{H.s.e) 



2e(l-2s^) (d-Lu+H 



X 



2swCw "V /g2=M2 



(118) 



where T^H+H-{,(f) is the sum of various bubble terms. 

These terms contain some combinations of Passarino-Veltman functions, such as: 

SE^{XY,q^) = AB22{XY,q^) + q^\B^{XY,q^) + B2i{XY,q^)-2Br{XY,q^)], (119) 



SEt{XY,q^) = 4S22(Xr,g2) + g2 B^{XY,q^) ^ B^^iXY.q^) 



(120) 



Here, we consider only the contributions which depend on q^, because r-^'"''''(if+iJ ) and 
p/m,z^^+^--) depend on the derivate of ^H+H-{(f)- Four types of bubbles are taken into 
account when calculating this "reduced" self-energy, which we refer to as f^H+H-iff)'- 

^h+h-{q') = tH+H-{VS,q') + tH+H-{SS',q') + tH+H-{ff:Q') + ^H+H-{xx':q')- (121) 
The VS bubbles contribute to TiH+H-{(f) with: 

tH.H-{VS,q') {e'SEf{H^, q') + ^'^^ ' ^l^^' SEfiHZ, q') 

4S|^ 4st, 



-SEt{h''W,q^) 



■>w 



4s^ 



SE^{A''W,q^)\. 



:i22) 
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The SS' bubbles contribute to T.H+H'il'^) with: 

light iff') ij=h^ 

+ gloHHBoiHH', q') + glo^HBoinh', q') 

+ gloGHMGH^ q') + gloGHBoiGh', q') 

+ g%GHBo{GA^,q')}. (123) 

The fermion and gaugino bubbles contribute to ^H+H-{q^) with respectively: 

+2M)il^,5o(//',g')} (124) 



and 



(126) 
obtained 



M^o + 5o(x.x;, g')}. (125) 

After having computed the full expression for T^H+H-{q^), one simply needs to calculate its 
derivative at q'^ — M% and then insert it into equations (117) and (118). 

5.3 Neutral Higgs sector 

The amplitudes a{"^ corresponding to final vertices with H^A^ or h'^A^ are: 

tTf^^'HH'Ayh'A^) rigi^^n iVf^^^^jH^A'/h^A^ 
a^^j,{H A/hA)- - ^— — X - . 

p/m,7(^o^o/;^o^o^ and r^*"'^(i/M7/iM°) depend on the final state. They are obtain 
by summing the contributions of various triangles and of neutral Higgs self-energy terms. 

For the photon, one has: 

pM,7('^o^oy^o^O'^ _ r^o^o/^o^o(lch) + r^o^o//i0^o(2) — r^o^o/;j0^o(2,pmc/i) 

+ r^o^o//jo^o(3/) + r^o^o/^jo^o (xxx) + r^o^o//jOAo(4ch) 

+ r]^o^o//j0^o(6/) + r]^o^o//i0^o(6ch) + r]^o^o//i0^o(4-leg). (127) 

For the Z boson, one has: 

r/in,Z(^0^0/^0^0) ^ r|o^o/,cMo(lch)+r|o^o/,„^„(ln) 

(2) — TfjOAo/i-fiAa{2,pinch) + r|'o^o//j0^o(3/) + 

ixxx) + r|oAo/;,o^o(x°x°x°) + r|»^o/,,o^o(4n) 

[4ch) + rfo^o/;j0^o(5) + rfo^o//jOAo(6/) 
[6ch) + r|„^o/ftOAo(6n) + r|o^o/ftOAo(4-leg) + 
{H.s.e) + r^HOAo/hOAo{H.c.t). (128) 





J- H"A"/hOAO 


+ 


^ H^'>A^'>/h^'>A^'> 


+ 


J- H^A^/hfiAO 


+ 


J- HOAO/hOA" 


+ 


■pZ 

J- H"A"/h"A" 


+ 


pZ 
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5.3.1 Tril-type triangles 

The Tril-type triangles contribute to T^^'^''^ {H'^A'^ / h'^A^) with charged terms only: 

rWAo^o(lch) = ±^x^[Za,]xC,{HWH). (129) 
The Tril-type triangles contribute to T-f'^'^'^^H'^A^/h'^A'^) with the following charged terms: 
r|o^o/,o^o(lch) = ±^x^^j^^[Zab]xC,{HWH) (130) 
and with the following neutral terms: 

r|o^o(ln) = ""'ri^sT^ I - a) CiiA^'ZH^) + cos'iP - a) C^iG^'ZH^) 

047r S^rC^/ 

+ cos^(/3-a)Ci(yl°Z/i°) -cos2(/?-a)Ci(G'°Z/i°)}, (131) 

rfo^o(ln) = { - a) C,{A'ZH') - sm'{P - a) C,{G'ZH') 

+ cos\p - a)Ci{A°Zh°) + sm\(3 - a) Ci{G° Z h°)} . (132) 

5.3.2 Tri2-type triangles 

With — \ "'" ^ ^ the contribution of the Tri2-type triangles is: 

I cw I &w loi' y = z 

rWAo^o(2) = ,^^[Z,,\xC^{WEW). (133) 
However, one must also take into account the pinch term: 



(2,pmc/i) = --^{Z^\xB^{WW,q'). (134) 



5.3.3 TriS-type triangles 

Let V be either the photon or the Z boson. The Tri3-type triangles contribute to both 
p/m,7(//o^o/;^o^O) r-'^^"'^(//M7/iM0) with two different terms. 

The first term, with fermion triangles at the final vertex, is given by: 
rWAo^o(3/) = - 16^2,1^^ Y.^lM)yA9vLf-gvRf) [Cs{fff) - M}C'i{fff)] . (135) 
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In the previous equation, the term yj depends both on the fermion in the triangle and on 
the final state. More exphcitely, for {qu,qd) or doublets, one writes yf as follows: 



Vf = 



^sinacot/S cosatan/?^ 



sin/3 



COS/9 



for //M", 



(cos q; cot/9 sin a tan /3\ ^ ,n.n 
^ ■ 
sm p cos /9 y 



(136) 
(137) 



The second term corresponds to chargino triangles at the final vertex. 
For H^A° final states: 

ijk 



+ My My 



^ik ^A°kj^H°ji ^ ^ik ^A°kj^H°ji 



C'MiXiXk) 



^ik^'^AOkjCROji + ^Yk^(^A»kj'^HOji ^siXiXjXk)} 



X] {[^ki^^AOjk^HOij + ^ki^(^AOjk^HOij CsiXiXjXk) 
ijk 



^ki^'^AOjk^W^j + ^ki^'^AOjk'^HOij (^siXiXjXk) 



Xfc 



^ki^C^Ojk^^HOij + ^ki^CAOjk^^HOij C'^iXiXjXk)}^ (138) 



For /i°74° final states: 

^loAoixXX) = ^ E { 



^Yk^^AOkj^^hOji + ^Yk^^AOkj'^h"ji ^■iiXiXjXk) 



ijk 



^ik (^A^kj^h^ji + ^ik ^A^kj^hOji ^siXiXjXk) 

+^x,^x, [OYk'-c^kjChOji + OYk^'c^okj^oj,] C'MXjXk) 

+M^M^^ Ofi/'c';{okjC^Oji + OYk^c%f.jC^0ji C'r^iXiXjXk)} 



■g;^ { [^ki^^Ojk^Oij + ^ki^^AOjk^Oij ^siXiXjXk) 

ijk 



+M^M^^ O]^f c^0jkCh0ij + O^/'c'^Ojkcj^Oij C'^iXiXjXk) 



VL^R 



+%^x/. 



^kf^AOjk'^hO'ij + ^k'i^'^A^jk'^hOij ^siXiXjXk) 

^ki^c^ojkC^Oij + Ol/^c^j^cj^^oij C'^iXiXjXk)}- (139) 



In addition, neutralino triangles at the final vertex give no contribution to r-l''''"''"'{H^A^ / /i°y4°) 
but they enter into the expression of r-^^"'-^(if°>l7^°^°)- 
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For final states: 



OL R L , /-/lOfl L 
ik ^AOkj^HOji ^ ^ik ^AOkj^HOji 



ijk 



^ik''^%kj''^HOji + ^ik''^AOkj''^HOji C'siXiX^Xk) 



C3(x°x°x2) 



Olfn^^.A, + 0L^n5o,,n^o,,] C^'(X°X;X^)}- (140) 



For the h^A^ final states: 

rroAo(x°X°X°) = i E { [O'^n^.p^io,, + O°fe^nio,,nfo,j Calx^xjx^) 



167r2 



ijk 



Aj A J 
Ai Afc 



167r2 



^ik'f^AOkj^hOji + ^^k^^'^kj^hPji ^slXi'XjXfc) 

^ife'^AOfcj^^Ojj + O^^n%^.jnj^0ji C'^iXiX^Xk)} 

^fcf^AOjfc^feOij + ^t^'^%jk'^h°ij C3(XiXjXfc) 
'^ki ^AOjk^hOij + '^ki ^A0jk^hOij\ '^sKXiXjXk) 

Olfn\o,A, + 0^n^^Api}\ C'Mx'.xl)}- (141) 



ijk 



A,; An 



5.3.4 Tri4-type triangles 

Let V denote either the photon or the Z boson. The Tri4-type triangles contribute to both 
p/m,7(//o^o//^o^O) and Vf'''^^{H^A^/hPA^) with charged terms, which are given by: 



(4ch) = gvwG X {gwHAognOGH C^iGHW) + gA^GHQwHH^ Ca{WHG) } (142) 



and 



r^oAo(4ch) = -^gvwG X {gwHAoghOGnC^iGHW) + gAOGHgwHhoC^iWHG)}. (143) 

The coupling constants gwHHO, dwHho and gA^GH depend on the charge of the particles at 
the vertex. Note that, in the triangles considered here, W and G carry the same charge. 
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Tri4-type triangles also contribute to rf^"''^[H^A'^/h'^A^) with neutral terms: 
1 



i HO AO 



(4n) 



167r2 



{ 



9ZZH0 



+gzzho 



gHORORogzHOAO Ca{H^H^Z) + QhOHORogzhOAO C4{H^h^Z) 
ghOROHogZHOA" Ci{h^H^Z) + gROhOhogzhOAO Ci{h^h^Z) 
+gzZRO gROAOAogZROAO Ca{H^ Z) + gROAOGogZROQO CiiH^G^Z)^ 
+gzZhO ghOAOAogzROAoC4{h^A^Z) + ghOAOGogzROGoC4(hPG'^Z)^Y (1^4) 



i hOAO 



(4n) 



[gzzH" gh"H"H"gzH"A" Ca{H^H^Z) + gH"h"h"gzhOAO C/^{H^h^Z) 

+gzzhO gROhOhogzROAO Ci{h^H^Z) + gnOhOhogzhOAO Ci{h^h^Z) 

+gzZRO gROAOAogzhOAO Ci{H^A^Z) + gROAOGogzhOQO Ca{H^G^Z) 
+gzZhO ghOAOAogzhOAoC4(hPA'^Z) + ghOAOGogzhOGoC4(hPG°Zj^Y (145) 



5.3.5 Tri5-type triangles 

The Tri5-type triangles contribute only to r-^^"'-^(if°A°//i°A°). 
1 



ROAO 



(5) = 



167r2 
1 

T67r2 



gzZRogZROAO C5{ZZH^) + gzZhogzZRogzhOAO C^{ZZhP 
gzZRogzZhogZROAO C5{ZZH^) + gzzhogZhOAO C^{ZZh^) 



(146) 
(147) 



5.3.6 Tri6-type triangles 

Two Tri6-type triangles contribute to both Tf^''^^{H^A^/hPA^) and Tf^''^^{H^A^/h^A^). 

Let V denote either the photon or the Z boson, we focus on r^o^o//j0^o(6/) first. The 
coupling of A^ to sfermions is proportional to the corresponding fermion mass and it is thus 
negligible, except in the case of third generation squarks, with sfermion mixing. 

For H^A^ final states, one has: 

rWo(6/) = H {gVbJ,k9H"hbjgA0bfbJ^^{^i^j^k) + gvWbkSAObibj^ 

ijk=l,2 

X] [gvtitkgHOiiijgAOifyCeititjtk) + gviiikgAOuijgROtjik^eiiiijtk)]- (148) 



ijk=l,2 



For h'^A^ final states, one has: 

^loAo{6f) =^ {gVbiH9h0bib^gA0bfbk^Sih^k)+gvbilu9A0b^^^^ 



ijk=l,2 



87r2 



X) {gviiik9h0iiijgA0ijik^&{titjtk) + gviihdAOiiijghOijh^ei'titjtk)] ■ (i49) 



ijk=l,2 
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The 6ch triangles also contribute to both r^*"'T(i/M7/iM0) and r^*"'^(i/M7/i°^°): 



r^o^o(6ch) = ^MK_g^,^^[g^^^c,{HGH)-gvGGCe{GHG)], (150) 



r^o^o(6ch) = ^^ghOGH[gvHHCe{HGH)-gvGGCe{GHG)]. (151) 

Since neutral Higgs or Goldstone bosons do not couple to a photon, the 6n triangles 
contribute to Vf'"''^ {H°A^ /h^A'^) only. 

For H^A^ final states, one has: 

rf o^o(6n) = ^-^ [gzhOQO [ghOHOHogHOAOGoCeih^H^G'^) + gH°h°h°gh°A°G°CQ{h^hPG^) 

+gZHOG° [gHOHOHogHOAOG°Ce{H°H°G°) + gh°H°H°gh°A°G°CQ{H^hPG^) 
+gZhOAO [ghOHOHogHOAOAoCeih'^H'^A'^) + gfjOhOhoghOAOAoCeih'^h'^A'^) 
+9ZH0A0 [9HOHOHogH<^AOAoCQ{H'^H'^A'^) + ^f/jOj/oj/o^i/jO^o^oCe (-f/'°/l°A°)j | 

~87r^ 1^^''°^° [gHOGOGogh°AOG°CQ{G^G^h°) + gHOAOGoghOAOAoCeiG^A^h^) 

+gZHOGO [gHOGOGogHOAOGoCe(G^G° H'^) +gHOAOGogHOAOAoCQ(G°A^H°) 
+gZhOAO [gHOAOA09hOAOAoC(i{A^A^h^) + gijo AOcoghO A^GoCi^i^^ h^) 

+gzHOAO [gloAOAoCeiA'A'^H'^) + ^|o^oGoC6(^°G'°i/°)] }. (152) 
For h'^A^ final states, one has: 

rfo^o(6n) = [gzHOGO [ghOHOHogHOAOGoC6{H^H°G^) + gHOhOhoghOAOGoCelH^h^G^) 

+gZhOGO [gHOhOh°gHOAOGoC6(h°H°G°) + ghOhOhoghOAOGoCeih^hPG^) 
+gZHOAO [gh°HOHogHOAOAoCQ{H°H°A°) + gHOhOh°gh°AOAoCQ{H°h^A^) 

+9zh0A0 gHOhOhogHOAOAoCeih^H^A^) + ghOhOhoghOAOAoCG{h^h^A^) | 

~87r^ [ghOGOGogHOAOGoCeiG^G^H^) + ghOAOGogHOAOAoCe{G^A°H°) 

+gZhOGO [ghOGOGoghOAOGoC6{G'^G'^hP) +ghOAOGoghOAOAoC6{G'^A^hP) 

+9ZH0A0 ^ghOAOA09HOAOAoCe{A^A'^H^) + gHOAOG^ghOAOGoCeiA^G^H^) 



+gzh°AO 



gloAOAoC,{A'A^h') + gloAOGoC,{A'G'h')] }. (153) 
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5.3.7 4-leg diagrams 

For H^A^ final states, the 4-leg diagrams give the following contributions: 

sin(/3 — a) f ^ 



r^o^c (4-leg) 



167r2 



0„2 



{ [Bo{HW, Mjj,) - Bi{HW, M|o)] 
+ [Bo(HW, Ml) - Bi(HW, Ml)] } (154) 



r^o^o (4-leg) = (4-leg) + I 

Cw 



X 



167r2 V4sfl,cf^.^ 
{ [So(>l°Z, M|o) - Si(A°Z, M|o)] 
+ [So(i/°Z, Ml) - Si(i/'^Z, Ml)] }. 

For h^A^ final states, the 4-leg diagrams give the following contributions: 

cos{P — a) 



rZo^o (4-leg) = 



167r2 



rfo^o (4-leg) 



-fl^lr^ (4-leg) _cos(/?-a) 



167r2 



,44 



.3 



X 



5.3.8 Neutral Higgs self-energies 

Before estimating the self-energy terms, let us first define several useful expressions: 



vi = sin 2a sin 2/3 2~ cos 2a cos 2/5. 

= COS 2a sin 2/3 -|- ^ sin 2q; cos 2/3. 



SEl{XY,q') 
SEliff.q^) 

SEl{XY, q\ a, b, c, d) 

SEl{XY,q') 
SEl{XY,q') 
SE^{XY, a, b, c, d) 



2q'B,{XY, q') - A{M^) - {q' + MI)Bo{XY, q^), 
2MfBo{ff,< 
8\^{ad + bc) 



2M]Bo{ff,q') + A{Mf) + q'B,{ff,q'), 



q'B,{XY, q^) + A{MI) + M^B^^XY, q^] 
+ (ac + bd)MxMYBo{XY, q^)}, 
B,{XY,q^)-Bo{XY,q'), 
2B,{XY,q^) + Bo{XY,q^), 
-9>[{ad + bc)MxBi{XY, q^) 

+ {ac + bd)MY [Bo{XY, q^) + Bi{XY, q^)] }, 



(155) 



{ [Bo{HW, M^o) - B,{HW, M^o)] 
+ [Bo{HW, Ml) - Bi{HW, Ml)] } (156) 



{ [Bo{A°Z, Mlo) - Bi{A^Z, Ml. 
+ [Bo{H''Z, Ml) - Bi(i/°Z, Ml)] }. (157) 



(158) 
(159) 

(160) 
(161) 

(162) 

(163) 
(164) 

(165) 



27 



In the following, all neutral Higgs self-energies and the Higgs tadpoles Tno/ho are 

written as the sum of various contributions coming from the gauge and Higgs sectors, fermion 
pairs, gaugino pairs and sfermion pairs (where we consider separately the unmixed case and 
the third generation squarks with mixing): 

= S(g+H) + S(//) + S(xx) + S(xV) + S(/7), (166) 
= r^o/;,o(g+H) + Tn^/hoUf) + Tho,Axx) + Tho/„o(x°X°) + Tno/AFf)- (167) 

a) self-energies: 



The contribution of the gauge and Higgs sectors is: 

1 



Si/o^^o(g+H) 



167r2 



+2glHHoSE^,{HW, g2) + 2gl^HoSE^,{WW, q') 



+2g^^H0 2Bo{WW,q' 



9ZZH0 



2Bo{ZZ, q' 



+gjjOHHBo{HH, g2) + g^.^^BoiWW, q^) + 2gjjoGHBo{WH, q^) 
+\ [glwBoih'h!', q') + gloHOHoBoiH'H', q') 

^ [9HOAOAoBo{A^A^,q'^) + gHQQOQoBQ{ZZ,q' 



+ 



+9hmoHoBo{H^h^, q^) + g]jOAOG°Bo{A^ Z, q^ 
e2M^cos2(/3-a) 



9„2 



B^{WW,q^) + ^B^{ZZ,q^) 



+ - 



\2A{M^) - 



+ 



2r2 L 



\2A{Ml) - Ml 



^2 



[(3sin2 2a - l)^(M^o) + Scos^ 2aA{M^o) 



cos 2/5 cos 2a A{Ml) - A{Ml] 



Ao2 



;i - v,)A{M^) + {1 + v,)A{M^h)] }• 



The contribution of the fermion pairs is: 

1 



^HOHo{ff) 



47r2 



E^/X(4o/)^X'^^2(//,0- 



/ 



The contributions of the gaugino pairs are 

1 



647r2 



XI ^BliXiXji ff, C-HOjv (^H°iji (^H°iji ^HOji) 



1 



1287r2 



I] '5E°(x°x?, n^oji, n^oij, n^oij, n^o^J 



(168) 



(169) 



;i7o) 
;i7i) 
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The contribution of sfermion pairs consists of two terms: 

1 



splight 



iff) 



167r2 



/ 



9hOhoUl + dHOROf^f^] MM})}, (172) 



SSoTo(//) 



167r2 



167r2 



167r2 



y=l,2 



E {4 + fe°HO/,/,^(M? )] }. (173) 



f=t,b 



b) self-energies: 



The contribution of the gauge and Higgs sectors is: 

E,o,o(g+H) = {gl^oAoSE',{A'Z, q') + gl^oGoSE',{ZZ, q') 



+2g'^HhoSE',{HW, q') + 2gl,ahoSE',{WW, q') 
+2g'^wf^o [2BoiWW, q') - l] + gl^ho [^BoiZZ, q') - 1 

+gloHHBo{HH, q') + gloccBoiWW, q') + 2gloaHBo{WH, q') 



+ 2 
1 

+ 2 



9hOHOHoBo{H^H^, q^) + gffit^offiBQ{h^h^ , q 



0!,0 2^ 



91pa°a°Bq{A^A^, q^) + gloQOQoBo{ZZ, q^) 

+9W^hoBo{H%\ q^) + ^?2„^„^„5o(AOZ, q^) 



0„2 



B^{WW,q^) + -^B^{ZZ,q^) 



+^(Mm^)-m^ 



s2 \ 



+ 



2c2 



2yl(M, 



Ml 



(3sin2 2a- 1)A(M|o) + 3cos=^2ayl(M^o) 
cos 2/? cos 2a U (Mj) - ^(M|) 



4c2 



(l + T;i)^(M^) + (l-^;i)^(Mi)] }. 



The contribution of the fermion pairs is: 

1 



W(//) = -^E^/x(4o/)'x5SO(//,a- 



(174) 



(175) 
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The contributions of the gaugino pairs are: 

1 



647r2 ^ 



1287r2 ^ 



E -^^s (X°X;, ^^0,,, n^o,,, nj:o,j, n^o,,). 



(176) 
(177) 



The contribution of sfermion pairs consists of two terms: 

sSo(/7)nght -T^E^i { H^LL + 9loM.] ^o(/7, .^) 



- [dhOHOf.f, + A{M})}, (178) 



■<r-theavy 



iff) 



167r2 



E {^ftOt,*, ^0 (tit,- , g') + fi/^Oft-b, ^0 (^i^j , g') } 



ij=l,2 



167r2 



167r2 



E {4 [9kOHoU,MMl)+9,o,of^j,A{Ml)] } 



E {4 WoHOf^fnMMl) + ^/.o,o/,/,A(M| )] }. (179) 



f=t,b 



c) Mixed if^/i'^ self-energies: 



The contribution of the gauge and Higgs sectors is: 

SHO/io(g+H) = ^gzH"A"9zh"A"SE^{A^Z,q'^) + gzH"G"9zh"G"SE^{ZZ,q^ 



+2gwHHogwmfiSE^,{HW, q") + 2gwGH^gwGh^SEl{WW, 
+2gwwHogwwho [2Bo{WW,q'^) - 1 
+gzzHogzzh° 2Bo{ZZ,q'^) - 1 

+gHOHHghOHHBo{HH, q^) + gHOGGgh^GGBoiWW, q^) 
+2gHOGHghOGHBo{WH, q^) 



+2 

1 

+2 



gHOHOHoghOHOHoBo{H^H^, (f) + gHOh°h°gh°h°h°BQ{h^h^ , q^) 



Oi,0 2^ 



gH°AOAogh°AOAoBo{A^A^, q^) + gHOGOGoghOGOGoBo{ZZ, q^] 
+gh°HOHogHOhOhoBo{H^h^, q^) + gHOA°G°gh°A°G°BQ{A^Z, q^) 

e^M^ sin(/3 — a) cos(/? — a) 



+ 



3e^ sin 2a cos 2a; 

8^2 2 



Bo{WW,q') + -^Bo{ZZ,q' 



A{Ml.) - A{Ml,) 



cos 2/3 sin 2q; [A(M^) - A(M|) 



"4.2 



A(M|)-A(M^)]}. 
30 



(180) 



The contribution of the fermion pairs is: 

1 



47^2 



J2Ni xc^HOfcj:ofXSE',{ff,q'). 



;i81) 



/ 



The contributions of the gaugino pairs are 

1 



2 ^ ^^siXiXj: C^Ojj, C^Oj^, C^Oj^, C^O^j), 

E -^^slXiX?, r^kj, n^HOij, n^noji)- 

ij 

The contribution of the sfermion pairs consists of two terms: 



(182) 
(183) 



ir^heavy 



iff) 



167r2 



167r2 



167r2 



E {^i/otit,- 9h0ui^ Bo {iiij , g^) + ^^os,s . QhObSj Bq (h,ij , g^) } 



(184) 



ij=l,2 



E {4 [^f.o,o/,/,A(M? ) +^^o,0/,/,^(M| )] } 



f=t,b 



d) self-energies: 



The contribution of the gauge and Higgs sectors is: 

E^o^o(g+H) {gl^^oAoSEl{h'Z,q^)+glH,^oSEl{H'Z,q^ 



+2g'^HAoSE',{HW, q') + 2g\oGHBo{WH, q') 
+gl.AOA^.B,{A'H\ q') + glo^o^oBoiA'h', q') 
+9hoaogoBo{ZH°, q^) + gloAOGoBo{Zh° , q^) 

+^ ([2A{Ml) - M^] + -V [2A(M|) - M 

cos^ 2(3 , , , ^2 N 



4c2 „2 



4s2 



1 + sin^ 2/? - ^ cos^ 2(3 



--w 



8s2 r2 



8s2 r2 



(3sin2 2/3 - 1)^(M|) + 3cos2 2/3^(Mj) 
cos 2/3 COS 2a [^(M^o) - A(M^o)] }. 



(186) 
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The contribution of the fermion pairs is: 

1 

64^ 



^A^Aoiff) =^ J2NlSEl{ff,q',c'Xof,-c'^of,c%,f,^c^AOf)- (187) 



The contributions of the gaugino pairs are 

1 

647r2 



Saoao(xV) = Y2L^E'^^3(x°X,^g^r^io,.,-r^k,■,r^io,,■,-^^^^^ (189) 

ij 



The contribution of the sfermion pairs consists of two terms: 

SS*.(/7) = -T7LEA'/{[w/./.+Sam«/,/J^(A^/)}. (190) 



16.^ ^ 



2 



16^' .,=1,2 



^ E {4 [^^o^o/,/,^(M? ) + 9aoaoM,A{MI)] } 



167r2 

f=t,b 



^ E {4 [w/./«^(^/,) + W/./.^(M? )] }. (191) 



167r2 

f=t,b 



e) Mixed self-energies: 

The contribution of the gauge and Higgs sectors is: 

S^oz(g+H) = _ f, , {Mzcos{P-a)sm{f3-a)SE',{H''Z,q' 



-Mz cos{p - a) sm{p - a)SEl{H'Z, q') 

+^Mz cos 2/3 cos(/3 + a) sin(/3 - a)SE^{A°H^, q^) 

-\mz cos 2/5sin(/5 + a) cos(/3 - a)5E°(A°/i°, g^) 

--Mzsin2/3cos(/3 + a) cos(/3 - a)SE°{ZH°, q^) 

+ -Mz sin2/3sin(/3 + a) sin(/3 - a)5£;^(Z/i°, g^)! (192) 
2 

The contribution of the fermion pairs is: 

W//) =^ E^/'5^6(//,?',^/ + «/,^/-«/,cio^,-cio^)- (193) 
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Here, Vf and a/ are defined as follows: 



1 



2s\yC\y 
1 



+ 



Vf = 



2 3'^. 

1 2 

2 + 3'^' 



a/ = + 



if / = g„ or 



AswCw 

Qf^-—^ — if/ = gd0r^. 

AswCw 



The contributions of the gaugino pairs are: 

Saoz(xx) = 8^ E {Mx, {c^AOjiO^f + c^o,,0,f ) [5o(x,X,-, q') + i?i(XiX,-, 



(194) 
(195) 



(196) 



Saoz(xT) 



1 



16,. E + <o,.Of ) [i^o(x°X?, + i^i(x^X;, 



Note that there is no contribution from sfermion pairs to S^02(g^). 
f) tadpole: 

The contribution of the gauge and Higgs sectors is: 

THo{g+R) {gHOHHA{M]i) + gnOGoAiMi 



w) 



+ 7;9hOhOhoA{M'}jo) + -gHOhOhoA{M^Q) 



2^ 



+ ^gHOAOAoA(M^o) + -gHOG°GoA(Mz) 



4A(M, 



2 ^ 2iV/? 



1 

■^gzzHO 



^eMw cos(/9 — a) 

Sw 

The contribution of the fermion pairs is: 



A{M^) + ^A{Ml) 



}■ 



The contributions of the gaugino pairs are: 

Tho{xx) = 



1 



^4o,,il^A(M|p, 

i 

8,2E^k.M,oA(M|o). 



47r2 
1 



(197) 



4A{Ml) - 2M| 



(198) 



(199) 



(200) 
(201) 
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The contribution of the sfermion pairs is the sum of two terms: 



167r2 



Y: {9HOuuA{MQ+9HouAiMO}- 



1=1,2 



g) Z;*^' ta(li)()lc: 



The contribution of the gauge and Higgs sectors is: 

T;,o(g+H) =--L {g^,o^HA{Mjj) + g^,OGGA{M^) 



1 1 

1 1 

+ 2^ftMOAoA(Mlo) + -ghOGOGoA{Ml) 

-gwwho [4A(M^) - 2M^] - ^qzzho [4A(M|) - 2M; 
eMw sin(/3 — a 



The contribution of the fermion pairs is: 



Tnoiff) 



47r2 



The contributions of the gaugino pairs are: 

i 



The contribution of the sfermion pairs is the sum of two terms: 



- YL^T.Nl{[9,ou.+9,of,j,]A{M})}, 



1=1,2 

h) Expressions of the renormahzed self-energies 



(202) 
(203) 



(204) 



(205) 



(206) 
(207) 

(208) 
(209) 



When calculating the effective contribution of the neutral Higgs self-energies to the pair 
production cross section at the one loop level, we must consider the renormahzed self-energy 
terms E, obtained by adding various counter terms to the unnormalized self-energies E. 
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Let us first define the Higgs field renormalization constants: 



5Zhi = 



5Zh2 = 



+ 



Mz 
tan/9 



(210) 
(211) 



They are used in the calculation of the various mass counter terms, together with the Higgs 
tadpoles and the parameter 5M^o^o defined as: 



5m1o^o = s^o^o(m1)-m2 



A" A" 



V dq^ J^,^^,^ 

Indeed, one has the following expressions for the mass counter terms: 

SMjjOHO = sin2(/3-a)5Mlo^o + cos2(/3 + a)Ezz(M|) 

- ""T^^J"^ ! + - «)1 '^Ho - cos(/3 - a) sin(/3 - Q;)r^o} 

+ M| cos(/5 + a) cos(/? - a) [6Zh, - SZh^] 
+ Ml cos^{P + a) sin^ P + (5Zh2 cos^ P 

SM^OhO = cos2(/3-a)5M|o^o + sin2(/? + Q;)Ezz(M|) 

+ ""T^^J"^ ! sm(/3 - a) cos(/3 - a;)rffo - [l + cos2(/3 - a)] T^o} 

— M| sin(/9 + a) sin(/? — a) — 6ZH2] 
+ Ml sm^{(3 + a) \6Zh^ sin^ /3 + dZn^ cos^ /9 



(212) 



(213) 



(214) 



5M|o^o = -sin(/?-a)cos(/3-a)5M^o^o -cos(/? + a)sin(/? + Q;)E2^(M|) 



-{ sin2(/3 - Q;)rffo + cos^(/3 - a)Tho] 



2swMw 
Ml sin a cos a [5^^/^ — SZh2\ 

Ml cos{P + a) sin(/3 + a) ISZh^ sin^ /3 + t^Z/f^ cos^ /3 



(215) 



There is no contribution from the neutral Higgs self-energies to r^^"'''^{H^A^/h^A^). As 
for their contributions to T-l^'""^ {H^A^ /h^A°), they can be derived using the renormalized 
self-energies, which are expressed as follows: 



^HOHo{q^) 

f 

Si?o/io(g^) 



^HORO {q^) + q^ SZffi cos^ a + SZh^ sin^ a — SM^ofjo , 
E/jO/jo(g^) -I- q^ SZffi sin^ a + SZh^ cos^ a — (5M^o^o, 
Sj?o/io(?^) + sin a cos q; [5Zh^ — 5Zh^] — SMj^of^o. 



(216) 
(217) 

(218) 
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For H'^A^ final states, one has: 
i HOAo[n.s.e) — 

SwCw 

For h'^A^ final states, one has: 

e cos(/3 — a) 



^1 n ■,'^HOhoiM]jo) 1 fdT,HOHo\ 



hO 



(219) 



SwCw 



tan(/3 — a) 



hOhO 



Ml, - M|o 2 



(220) 



5.3.9 Neutral Higgs counter terms 

For H^A^ final states, the neutral Higgs counter terms give the following contribution: 

e cos(/3 — a) 



{H.C.t) 



2s\yC\y 
esin(/3 — a) 



[cos /3 sin /3 + cos a sin a] {5Zh2 — SZh^ ) 
^cos^ a + sin^ /9)5Zhi + (sin^ a + cos^ /3)5^i/2l ■ (221) 



For h A final states, the neutral Higgs counter terms give the following contribution: 



6 sinf /3 — cx] 

(H.C.t) — [cos /3 sin /3 — cos a sin a] {5Zh^ — 5Zhi) 



+ 



2sy/Cw 

e cos(/3 — a) 
2s\yC\y 



(sin^ a + sin^ (3)SZh, + (cos^ a + cos^ P)SZh2 . (222) 



6 Contribution of box diagrams 
6.1 Diagram structures for boxes 

Several useful expressions are needed when estimating the contributions of the box diagrams. 
The particles inside the box are ordered clockwise and have internal masses mi, m2, m^, 
1714 starting with mi between the and junctions. In the following, we will use the 
Mandelstam variables s, t and u, which can be expressed as a function of and of the 
masses Mi and M2 of the two outgoing Higgs bosons (i.e. H^H" or H^A^ or hPA^): 



s — q , 
1 

t 



+ 



^(Mi + Mt-s) 
s cos 9 



2 \ 



^ M2 + Mi >| M2 + Ml 



M2 - Ml 



1 - 



M2- Ml 



(223) 



(224) 



1 



u = -{M^ + Mt-s) 



scosO 



2 ^^ 



1 + 



M2 + M1 



1 - 



M2 + Ml 



1 + 



M2 - M/ 



1 - 



M2 - Ml 



(225) 
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Note that u and t have an angular dependence and, as a result, so does the contribution of 
the diagram boxes to the one loop cross sections. 



Let i and £' be the momenta of the incoming electron and positron, respectively. Let Pfi 
and Pf2 be the momenta of the two outgoing Higgs bosons. When not otherwise specified, 
£', Pfi, Pf2 and £ are oriented clockwise around the box. Then, one has: 



p2 
P/lP/2 



Ml 
Ml, 

s - {Ml + M|) 



t-Ml 



u-Ml 



(226) 
(227) 

(228) 
(229) 
(230) 



a) Box7 structures: 



+ -P/l(-D222 — -D223; 

0^22) - 2£'Pf2iDn3 



Pf 2(^^332 — -D333 
-D^,2)-2£'Pf^(D 



+ 



6(-Doo2 — -D003 

2PflPf2iDs23 — D322) — 2i''P/2(-Dl33 ~ -D132) — 2£'Pfi{Di23 — D122 
4 [P|iD22 + P/2^23 + 2f P/1D24 + 2£'Pf2D25 + 2PfiPf2D2Q + 4^27 

4 \P]2Dn - Pj.Du - 2£'Pf,Dn 



(231) 



b) Box8 structures: 



— 6(1^002 — -D003) + P/l(-D222 — -D223) + P/2(P'332 — D'SSs) 

- 2Pf^Pf2{D^23 - D322) - 2£'P/2(/^133 - ^132) - 2f P/i(Di23 " ^^122) 

+ (2f P;i + P}^)D22 - (2f P/2 + 2P^iP^2 + P;2)^23 " (2f P/2 + 2£'Pf,)D2, 

+ {2£'Pf2-2£'Pf,-2P],)D26-2D27, (232) 



D 



12 



D 



13, 



c) Box9 structures: 



d) BoxlO structures: 



D" = Pfo + P'i2 - P'la- 



Vg = Pi2-Pl3. 



Pio = P'12 + P'o, 



P'lo = ^12- 



(233) 
(234) 

(235) 

(236) 
(237) 



Finally, the crossed functions T>j are obtained from the Vj functions by making the following 
changes: P/i P/2, i.e. t ^ u and M^ M|. 
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6.2 Charged Higgs sector 

For e+e" H~^H~, the box diagrams give the following one loop contribution: 



2e2 



-rA'£${H^H-)-rA'£$\H^H-)]. 



(238) 



6.2.1 Box7 diagrams 



The amplitude A'f^^{H~^H ) is obtained by summing various contributions with gauge and 
Higgs bosons inside the box: 



A'rS{H+H-) 



a: 



+ 
+ 



a. 



8s4 



sm\p - a)PLV7{uWH^W) 
cos^{(3 - a)PLV7{iyWh'^W) 



PLVr{iyWA°W) 



2swCw , 



(I - 2s 



w 



\ 2sy^rCw 



X 



X 



9lPl + QrPr 



9lPl + GrPr 



+ 



(I - 2s 



2 ' 
W 



\ 2s\YC\y 



X 



2s\yC\y 
gjPL + qIPr 



X al^{Vr{e^HZ) -V,{ejHZ)} 
X al^{Vr{eZH^) -V,{eZH^)} 
X al^{Dr{eZHZ) -VrieZHZ)} 



+ [Pl + Pr] X al^{Vj{e^Hj) - V,{ejH^)}. 



(239) 



6.2.2 Box8 diagrams 



The amplitude A^£^^{H^H ) is obtained by summing two types of box diagrams, which 
have gauginos inside the box, I'xx'^x ex°xx°: 



L,R 



A'^%'h- ii>Xx'x) + ATA- {ex' XX 



boxS /'~~0~~0> 



(240) 



For the C'XX^X boxes, since there is no right-handed sneutrino in the MSSM, only a left- 
handed term is considered: 



ATh-{^lXX"x) 



167r2^ 



Y: Zt*Zt,X {M^M^^C^H\AkPli^LX^x]Xk) 



W ijk 



+M^M^ocf;,A,^V',{uLX^X%) 
+4*A^s{i>LM%)}- (241) 
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For the ex^XX^ boxes, one has both left-handed and right-handed terms. The left-handed 
term is given by: 

+M^^M^o4*,4^.,P^(eiX-X,-X°) 

+4;^4,fc^8(eLX-X,Xfc)} 

Y^{Z^i*sw + Z^*cw){Z^kSw + Z^^cw) X 



+M^„ M^^, c^,,4^.fcP; {eLx'lx.xl) 
%° 4}i4,-fc^8 (e-LX ■ X,-Xfe) 
+4)iC%k'D,{eLX%xl)} (242) 



while the right-handed term is: 



A'SfH-{eax'xx') =^jVE^iT^itx {M^oM^o4;,4,.,P^'(e~^x°X,X°J 



/ 



%^x/H;.4,-;t^8(eiiX-X,X^) 
+M^.M^ocf;,4,.,P^(e,^x'x,X^) 

+4}.4,fc^8(ei?X-X,X°)} 



8^2,2^ E ^iT X { A^x°^x^4 -.^..^s (ei^X°X,X°.) 



+M^oM^^4;,4^.,D;(g^x-x,Xfc) 

+^x,^x2ci;-«4,fc^^8(ei?X°X,Xfc) 
+4;i4,fe^8(ei?x° X.Xfe) }■ (243) 



6.2.3 Box9 diagrams 



There is no right-handed contribution from box9 diagrams. As for the left-handed amplitude 
A^£'-^^(H'^H~), it is obtained as follows: 



e2 



e2 



E l^il'^k.,^9(x.^LeL^L). (244) 



39 



6.2.4 Twisted boxlO diagrams 

Twisted boxlO diagrams contribute only with a left-handed term: 

Y,zUz^;sw + z^;cw)xgH, 



X 



[^x? 4iPw{eLX%yL) + M^^ 4,,-:Pio(eLX°X.^>L)] ■ (245) 



6.3 Neutral Higgs sector 

For e^e~ — > H^A^ /h^A^^ there is no box9 diagram, so one has: 

^60.(^0^0/^0^0^) {A'>£%{H^A^/hPA'')^A^£${H''A^lhPA'') 



^^60x10(^0^0/^0^0) I (246) 



6.3.1 Box7 diagrams 



In the neutral Higgs sector, there is no right-handed contribution from box? diagrams. The 
left-handed amphtude A^f'^^H^ A^ / hP A^) is obtained as follows: 

^60.7(^0^0/^0^0) = ^ X [Z,,] X {v^{vWHW)+V^{vWHW)]. (247) 



6.3.2 Box8 diagrams 

Both vxXX ^iid ^x°x°x° box diagrams contribute to the amplitude A^£^^[H^ A^ / hP A9): 



^6o.«(^u^u/^o^o) ^ P^,«[A^%lo/,o^o(^^XXX)+^»/.o^o(exW)J- (248) 
The I'XXX boxes contribute only with left-handed terms: 

A'A i^LXXX) = tSt^ E ZtCZtk X {M^M^^c%^,c^HOkPl{^LmjXk) 

+M^M^^c%^^4okPsiHx,XjXk) 



..^2,2 E^it^ix {A%M,^4o,,ciofc,P^'(z>iX,XiX.) 



+Mx.M^,c^0jiC%kP'8{uLXiXjXk) 
+CH0jiC%kP8 {i>LXiXjXk) } , (249) 
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W ijk 



+^X, ^X, CAOjiChOkjT^'sii^LXiXjXk) 

+^x, ^x, CAOjiChOkj'D's {i>LXiXjXk) 
+CA°jiCh°kj'^s{i>LXiXjXk) } 



e2 



+MiM^.C^OjiC%f,p'^{i>LXiXjXk) 

+M^M^,ChOjiC^OkjV'siOLXiXjXk) 
+Ch0jiC%kj'^8{i>LXiXjXk)}- (250) 
The ex°x°x° boxes contribute with both left-handed and right-handed terms: 



{^x?^x'^"'A",;"4o,fcP;'(eiX-X;X°^ 
+M^oM^on:^o,,ngo,feP^(eiX?X?Xfe) 

g2 

+ Z2i*Cw){Zi^Sw + Z2kCw) X 

{M^oM^.mgo,,nio,,:P^'(gLX-X°X°) 
+M^oM^on|o,,n^o,,:P;(gLX-x;Xfc) 
+M^oM^on^o,,nio,feP^(eiX°X?Xfe) 
+^io,,.<o,.,P8(eLX-X°X°)} 



--5-^P«^ X {M^oM^onio,,.n|o,.,D;'(e«x°x;x°J 

"^w ijk 

+M^oM^onio,,n^o,.,P;(e«X-x;Xfc) 
+M^oM^.m^„,^.ng„^.,I)^(g«X-X;Xfc) 
+n:^o,,nio,-,I)8(ei?X°X;Xfc)} 



--^P^ Y: ZlZl* X {M^oM^onio,,.n5o,.,P^'(e«x°x;xD 

ijk 

+M^oM^on^„,^.n^o,,P^(gKX°X;Xfc) 
+M^oM^on§o,,n:^o,feP^(eflX°X;Xfe) 



+^fo,,.n^o,.,P8(eiiX-X;X^)}, (251) 
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PlY.{Z^:sw + Z^*cw){Z^uSw + Z^cw) X 



X 



ijk 



{M^oM^on^o,,4„,fcI?;'(gLX-x;Xfc) 
+M^oM^on5o,,nfo,,P^(g^X-X?X°) 
+M^oM^on:^o,,n^o,.,P;(ei,x°X°Xfe) 
Moi,<o,fcP8(eLX-x;x°)} 

e2 Af* AT JV 

{M^oM^„nfo,,nio,fc:D^'(gLX-X?X°) 
+M^oM^onfo,,n5o,fcP;(gLX-X?Xfc) 
+M^„M^„n^o,,nio,fcP;(giX-x;Xfe) 
+i^mji^A0jk'^8 (eLXi X?Xfe) } 

°^ iik 

+M^oM^onio,,n^o,feX);(eflX°X;X°) 
+M^oM^on2o,,n^o,.,P^(e«X-x;x°) 

{M^oM^on^o,,.n5o,.,P^'(e^X°X;X^) 

+M^oM^o4o,,nio,fc:P^(gfiX-x;x°) 
+M^oM^on^o,,n5o,,P^(e~^,X°x;Xfe) 
+^ft0ij^A0jfe^8(eRX°X?Xfe) } ■ 

6.3.3 Twisted boxlO diagrams 

Two types of twisted boxlO diagrams must be considered: i^xX^ ex°x°e. 
The I'XX^ boxes contribute with left-handed terms only: 

^2 



N ryN* 

Ik 



X 



ijk 



CAOji'^w{i>LXiXji>L) + C^OjP'ioi'^LXiXjj^L) 



VL^L 



X 



CA0ji1^10{j>LXiXj^L) + M^. C^jP^Q^yLXiXj^L) 
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The ex°x*^e boxes contribute with both left-handed and right-handed terms. Writing all 
these terms into one single expression leads to: 

„2 



ij 



A»(ex°x'^e) 



[m^o nio,,Pio(eflX°X;eii) + ^AO.i^^ioleflX^xJeii)] , (255) 

Pl J2iZu^w + Z^*cw){Z^^sw + Z^.cw) x ^,0,,,, x 
[m^o n:^o,-,Pio(eLX°X?eL) + M^o n:^o,-,Pio(eLX°X?eL) 



X 



[m^o nio,,Pio(eKX-X;eK) + M^o n^o,,Pio(eiiX-X;efl)] . (256) 



7 Conclusion and outlooks 

In this paper, we discussed all electroweak one loop contributions to the pair production cross 
section for charged and neutral Higgs bosons in e^e~ collisions, in the theoretical framework 
of the MSSM. The one loop amphtudes of initial vertices and self-energy, of 7 and Z boson 
self-energies, of the corresponding counter terms, of final vertices and Higgs self-energies and 
of box diagrams are respectively given by equations (17), (82), (83), (98) and (238) for the 
charged Higgs sector, and by equations (18), (85), (86), (126) and (246) for the neutral 
Higgs sector. The left-handed and right-handed amplitudes of all these electroweak one loop 
contributions are: 

• in the charged Higgs sector: 

+ af:^{H^H-) + a%^{H+H-) 

+ a%{H+H-), (257) 

• in the neutral Higgs sector: 

ai;7(i/°^7/i°^°) = a™«(/J°A7/?°A°) 

+ af%{H^A^/h''A'^) + ag^(i/°^7/i°^°) 

+ aS(/J°A7/iM°) 

+ a^°^(ii'M7/i°A°). (258) 



43 



The differential production cross section for charged or neutral Higgs bosons at the one loop 
level can then be calculated as follows: 



do ^ 7r«|^ 2 Q-s ^ [|a£orn|2 ^ 2|af °""ai'°°^| + |a^°""|' + 2|a|°™a],'°°^|l . (259) 

dcOS^ 8g2 ^ ' V ^ \ \ L L \ \ R \ \ R R \\ \ ) 

In the previous equation, af'^" is the Born amplitude of equation (8) or (9). As for /S^, 
it stands for the velocity of the Higgs bosons, see equation (12). After integration over 
cos 6* (which appears in the contributions of the box diagrams), one obtains the total pair 
production cross section at the one loop level. Note that, in the case of the tree level cross 
sections for H'^H~ and H^A^ + hPA^, there is a direct dependence on Ma only and not 
on the other MSSM parameters. However, after having taken into account all electroweak 
one loop contributions, this is not true anymore. Indeed, Oj^"^^ depends on other MSSM 
parameters than just Ma- 

A C++ numerical code has been developed in order to calculate accurately all one loop 
electroweak contributions and, in turn, to compare the pair production cross sections for 
MSSM charged and neutral Higgs bosons at tree level and at the one loop level. The relevant 
Feynman diagrams are computed by calling the suitable functions in the Loop Tools 2.1 
library [16]. The input of the code, in standard notation, is the following: tan/3, /x. Ma (the 
mass of the Higgs boson), the gaugino parameters Mi and M2, the scalar mass scale M5, 
the sfermion mixing matrix parameters A„ and A^. A possible reference for these parameters 
is [15]. This input requires a preliminary pre-processing using the FeynHiggs 2.1 [17] code. 
A subset of these parameters is then fed into FeynHiggs in order to compute the masses of 
the Higgs bosons hP, H"^, as well as the mixing angle in the neutral sector a. These 
additional parameters do indeed appear in the analytical expressions described in the text. 
The output of the code is the cross section for the various processes under consideration. 
We have successfully checked that the MSSM Higgs bosons pair production cross section 
computed by our code at the one loop level remains stable against UV divergences, both 
in the charged and neutral Higgs sectors. Also, we have checked that the variation of the 
computed one loop cross section with (f agrees with our expectations. 

Note that, in this code, we have included all virtual contributions involving particles having 
electroweak interactions in the MSSM, but we did not treat pure QED effects, such as Initial 
State Radiation (ISR) and Final State Radiation (FSR) . The reason is that these effects may 
depend on the characteristics of the detectors (for instance, they need specific kinematical 
cTits) and some specific codes exist in order to treat them. Nevertheless, in order to be able to 
test electroweak symmetry properties at high energy, in particular those of supersymmetric 
nature, which is indeed the purpose of this work, we have included the virtual photon effects, 
but with a photon mass set to Mz in order to keep these effects finite. In order to have the 
complete (observable) contribution including QED effects, one should compute the following 
combination: Our contribution + ISR + FSR + virtual soft photon with zero mass — virtual 
soft photon with Mz ■ This combination should be calculated at the level of the codes that 
include the ISR and FSR effects. 
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Appendix A: Vertices and couplings 

Gauge - Fermion - Fermion 

Let Qf and T| be respectively the charge and the third isospin component of the fermion /, 
then the couphngs of gauge bosons to left-handed and right-handed fermions are: 

PLg-,ff = Qf and PjiQ^ff^Qf, 

PLgzff = — and PRgzff = -Qf — , 

PL9wff' = — ^ and PrQw}}' = 0. 

Note that, in this paper, we have also used a simplified notation for the couplings of 7 or Z 
to fermion pairs: 

gvLf = Pigvff, 

gvRf = Pngvff, 

where V stands for either 7 or Z. 
Gauge - Gaugino - Gaugino 

Ol'' = -eAj and Of = -e6,j, 



e 

Ol^ = — 



w) 



and 0,^." = 

2swCw 2swcw 

„(7N*yN yN*yN\ 7N yN* yN n'N*\ 

qOL ^ ^y^ii - ^OR ^ el^4i^4j - ^3i^3j ) 

'■^ 2swCw 2swCw ' 

~ ^ [ 72 ) ~ ^ [ ^ 7^ ) ■ 

Here, the Zij terms correspond to the various elements of the unitary mixing matrices of 
the charginos and neutralinos. They are derived from the diagonalization of the MSSM 
gaugino mass matrix, see for instance reference [18] for details. 

Gauge - Sfermion - Sfermion 

The couplings of gauge bosons to unmixed left-handed and right-handed sfermions are: 
3zhh- ^TT- gzf,f, = eQs — , 



SwCw ^JHJH 
^0 
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Let 6 J he the mixing angle of the sfermion / (generally a third generation squark). Let 
us also define = cosOj and = sinOj. The coupling between a gauge boson and two 
sfermions with mixing is then given by: 

+ cV- - 







9zfJ^ 




9zfJ^ 


= 9zfJ, = 




9whh ^ 




9whf2 ^ 




9whf2 ^ 



J'9wfLfi^^ 
'i'9wfLf'i^^ 
f^f'9wfj'^^ 
9whi[ ^ ~^f'^f9wfLf'i^- 

Gauge - Gauge - Higgs 

gzZHO = 2~ ~ ^) 9ZZhO = 2~ ^^^(/^ ~ '^)' 

gwww^ = cos(p - a) and Qwwho = sm[fj - a), 

Sw Sw 

g-iWG = eMw and gzwG = -eMw — - 

Cw 

Gauge - Higgs - Higgs 

6 6 

= sin(/? - a) and g^z/jOAo = +7. cos(/? - a), 

ZSwCw ZSwCw 
6 6 

= +7^ cos(/? - a) and g'zfeOGO = +7^ sm{(3 - a), 

ZswCw ZswCw 

6 6 

5'ty±if±ifO = + X - — sm{(3-a) and g'wiifift" = - x 7^ — cos(/3-q;), 

ZSw ZSw 

g g 

fl'H^iGiifo ^ -Qw X - — cos(/3 - a) and g'w^iGifeo = - x - — sin(/? - a), 

^Spi/ ZSw 

e 

gwHAo = +- — and gwHG" = 0, 
2sw 

g 

gwGAO = and gwGG° = +7^ — , 

ZSw 

g-yHH = -e and qzhh = ~e— , 

ZswCw 



g^GG = -e and qzgg = -e 



w 



2swCw 
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Higgs - Fermion - Fermion 



The coupling constant between a Higgs boson and a fermion pair is proportional to the 
mass of the fermion(s). Thus, only the third generation quarks are usually considered. 



In the charged Higgs sector: 



^/2swMw 



Mt cot/3 and c? 



b-^tH- 



Mfetan/?, 



w 



• for i ^ 6 H+: ci\,„+ = ^ ^ ^ ' ^ 



Mb tan P and cf_^^jj+ — 



In the neutral Higgs sector, the left-handed coupling constants are: 



Mt cot/3. 



• c 



• ^hH — 



eMt sin a . 

X — — - and c^ofe 



eM, 



2swMw sin (5 



h cos a 

X 



i cos a ^ 
X — — - and 



2swMw sin /? 



//Oft 



2swMw cos /3 ' 
eMft sin a 

+ - X 



2swMw cos /3 ' 



X cot/3 and c^o^ 



X tan /3. 



2swMw 2swMw 
As for the right-handed couplings constants, one simply has: 



• C 



Higgs - Gaugino - Gaugino 



In the charged Higgs sector, there are two types of vertex to consider: x° — > xt H~ and 
Xj' x^ . Let us choose the case where i and j label respectively a chargino and a 
neutralino, then the corresponding left-handed and right-handed coupling constants are: 



• c 



Hij 



esin (5 

SwCw 



1 



ecos/3 

SwCw 



For the other vertex (where i and j label respectively a neutralino and a chargino), one 
should instead use the left-handed and right-handed couphng constants c§*j and c^*j, 
respectively. 
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In the neutral Higgs sector, one must consider the couphng between a neutral Higgs boson 
and either two charginos or two neutralinos. 



\/2sw 

e 

\/2sw 

e 



For the coupling between a neutral Higgs boson and two charginos: 

6 r 1 

\cosaZ^^Ztj + s\naZ^^Z:^^\ and cfoij = c^o^j, 
- sin aZ-Zt- + cos aZ-Z+^ and c^o^^ = c^o^^, 
sin (3 Z^iZtj + CO& (3 Z{^Z:^^ and 0%^^ = -c^lji- 



V^Sw 

Let us now consider the coupling between a neutral Higgs boson and two neutralinos. For 
the left-handed components, one has: 

X { (cos aZ^. - sin aZ^-) {Z^^sw - Z^, cw) 

+ (cosq;Z^ - sin aZ^i){Zf^jSw - Z^jCw)}, 
{ (sin aZ^j + cos aZ^j) {Z^^ sw - Z^^ cw) 



■X 



L 

rijyOij 



■X 



-|-(sinQ;Z^ + coso;Z^)(Z^siy — Z|^cvi/)}, 
2^^^^ { (sin PZ^^ - cos pZ^^){Z^,sw - Z^, cw) 

+(sin/3Z3^ - cos (5Z^:){Z^^sw - Z^^cw)}. 
As for the right-handed components, one simply has: 



"'HOjii 
''^hOjii 



Higgs - Sfermion - Sfermion 

Let us first consider the light unmixed sfermions. Their coupling constant to the charged 
and neutral Higgs bosons is not proportional to the mass of the corresponding fermion(s) 
and it can thus not be neglected. 

In the charged Higgs sector, if / and /' represent respectively up-squarks and down-squarks 
of the first and second generations, or sneutrinos and charged sleptons, one has: 

eMw . . n 



In the neutral Higgs sector, one has: 



eMw „ 2 



{Tf - QfSw) cos{a + (5) and g^o 



SwCw 

9h0fj\^^^iTf-Qfsl,)sin{a + p) and g^ofj 



eMw „ 2 / , m 

SwCw 
eMw ^ 2 • / , r,^ 
SwCw 
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Note that the couphngs between and a pair of hght unmixed sfermions are proportional 
to the fermion mass and are thus neghgible. 

Let us now consider the heavy sfermions, i.e. the third generation squarks, and, in a first 
step, let us assume that there is no mixing. The coupling constants between the MSSM 
Higgs boson and a pair of unmixed heavy sfermions are given below. 



For the charged Higgs bosons H, one has: 



_ eMw 



9m 



Sw 



V2 



MHaji 13 + Mf coil3' 
sm2/5 * 



For the neutral Higgs boson A^, there are only off-diagonal terms: 

SaHJ^ = aVt^i^ = and g\,-^j^^ = g%~^j^ = 0, 

9 AH J, = 9 AH J, = [-/^ - Atcotp] , 

9A0iJ^ - 9aoI,1, - [-/^ - Atan/3] . 

For the neutral Higgs boson one has: 

eMw (I 2 ^\ eM^ sing 

eMw /2 2 \ . eM? sin a 



= [r-^ ) ^ - swMw sin (3^ 

Q Q eMt / — /i cos a + AfSinaX 

9hHJ, = 9H0i,i,--^^^Mw [ ^ j 



(1 eMiy /I 1 2 \ / ^\ sM? cos a 

"i^i sh/c^ V2 3 / swMw cos p 

eMvi/ /I 2 \ / ^\ ^^b cos a 

9h% 1 = 2- U*i^ cos(q; + /3) — -, 

"r^r swCw V3 / swMw COS p 

Q Q eMft /— /isina + AftCOsaX 

9h%J>^ = ^H0SA^~2sw^MvK i ^^^^ J' 
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For the neutral Higgs boson h^, one has: 



n eMw /I 2 2 \ . / n\ eM? cos a 

n eM^y /2 2 \ . . ^. eM? cos a 



Q Q eMt ( —fj, sin a — At cos a \ 

n eiV'/vi/ /l 1 2 \ • / ^\ eM^^ sin a 

Q eM^y /l 2 \ • / ,3\ sin a 



^'^"v. - -;;^lr^J^^"^"+^)+.^M^cosr 

Q Q eMb /— //coso; — AbSinaX 

9h0bJ^ = ^/.°6A"~2svfM^\ ^^^^ )■ 

Let us know take the sfermion mixing into account and let i?* and R'' be the rotation matrices 
for i and b squarks, respectively. If /{'g = then: 

^ = 4^with4 = (^ J/). 

In the charged Higgs sector, one has: 

' i'j' ' ' 

Similarly, in the neutral Higgs sector, for / stands for either i or b, then one obtains the 
following coupling constants: 



i'j' ' 



Note that, in the case of the neutral boson A^, the coupling to a pair of sfermions is the 
same with or without mixing. 
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Higgs - Higgs - Higgs 



eMw \ cos 2Pcos( (3 + a) ,^ , 
gn^HH = TT-^ cos(/3 - a) 

eM^ji/ [cos 2/9 sin (/3 + a) 

ghPHH = 772 ^ ~ ^) 

gn^GG = -7r~~^ ["^o^ 2/3 cos(/3 + a)] , 

^/lOGG = 2s^5~ 2/3sin(/3 + a)] , 

sin 2/3 cos(q; + /3) 

9 



.(/3-a) 



gHOGH = 7: 

2sw 

eMw \ . 
ghOGH = -7, cos(/3 - a) 

^ eMw 

gAOG±H± = X -7, ) 



r2 



sin 2Psm{a + (3) 



(-2 



— cos zcK cosyp + q; j , 

2swCtia 



ghOHOHO = — [2 sm 



g HO HO HO — — 7S 2~ 

S'/iO/iO/iO = 2-cos2Q;sin(/3 + a), 

ZswCw 

^^'^ [2 sin 2q; cos(/3 + a) + cos 2q; sin(/3 + a)] , 

IswCw 

gnOhOhP — ~ 7. '~2' [2 sin 2q; sin(/3 + a) — cos 2a cos(/3 + a)] . 

eMw / n \ 

gnOGOGO = -7^ 2~ 2/? cos(/3 + a), 

eMw 

ghOGOGO = — cos2/3sin(/3 + a), 

g HO AO AO = TT—^ cos2/3cos(/3 + a), 

ZSwCw 

eMw ■ m \ 

ghOAOAo = -7^ 2~ COS 2/3 sin(/3 + a), 

gnOAOGO = 7r~~^ 2/3 cos(/3 + a), 

S'ftOAOGO = - ^^^2 sin 2/3 sin (/3 + a). 

ZSwCw 



Higgs - Higgs - Sfermion - Sfermion 



For light unmixed sfermions, the couphng constants are not proportional to the mass of 
the fermion(s) and can not be neglected. The coupling constants for the heavy sfermions 
are then obtained by adding a term proportional to the mass of the corresponding fermion(s). 
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If / is a slepton {£ or or a squark from the first or second generation (g), one has: 



dnOHOfufii 



0„2 



cos 2a 
COS 2a 



2s^ 



w 

2 



W 
2 



9H0h0fRfR = ^ 



0„2 



2 



COS 2a 
COS 2a 
sin 2a 
sin 2a 



(Qf^w) 



(Qf^w) 



r2 



9s2 



9A0A0fRfR - r, 2 



COS 2/3 / 3 2 \ 

COS 2/? / 2 \ 
[QfSw) 



r2 



For the third generation squarks (stop and sbottom), one has: 

f Mt sin a 

9H0H0iL,RiL,R = 9hohOul^rul,r - TTT 



9h0hHL,Rh 



QhOhOu 



2slr \Mw sin p J ' 
/ Mt cos a \ ^ 



sin 2a / M. 



9H0H0i,,Rh,R - 9H0H0u,,RU,,R ~ sin ^ ) ' 



9A0A0iL,RtL,R - 9A0A0ui^^rUl,R " 77^ 



/ Mt COS p 



and 



2sf^ VMvKsin/?^ 



/ Mucosa 



9H0H0h^RH,R - 9H0H0d,,Rd,,R 2s'w[MwCOsP ' 



/ Mfesina 



- 9hOhoi,^Ri,^R 2sI,[MwcosP) ' 



e^sin2Q; / Mb 



9H0h%,^RH,R - 9HOhod,,Rd,,R 4^2^ [mw cos P J ' 



fl'AOA06i,fl6i,fl - 9A0A0dL,RdL,R 9^2 



/ Mftsin/? 



2s^, \Mw cos P 
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Appendix B: Passarino-Veltman functions 

The calculation of one loop Feynman diagrams can be performed by combining propagators 
using the following formula: 



/ dxi ■ ■ ■ dxnSixi^ hx„-l) 

io 



• • • A?r r(Q;i) • • • v{an) h " ' ' " ' (Aixi + • • • + ■ 

However, it is often convenient to reduce each one loop diagram to the sum of standard 
contributions, the so-called Passarino-Veltman (PV) functions. 

B.l Standard definitions 

Let us define the 1, 2, 3 and 4 point functions according to: 

d^k 1 



J I 

{Do,D^',D^'\D'"'P}{abcd) = 
where the denominators are: 



d^k{l,k'',k''k''} 
d^k{l,k^',k^'k''} 
d^A;{l,A;^,A;^r,A;^rA;''} 



Ni = k'^ -ml + ie, 

N2 = {k + pif - ml + ie, 

N3 = {k + pi + p2)^ - ml + ie, 

N4 = {k + pi + p2 + ps)"^ - ml + ie. 

Here, all integrals are kept D-dimensional. However, the rest of the calculations will be 
performed in four dimensions. 

In one loop diagrams, the following conventions are used: 

• the external momenta pi„.N are oriented clockwise, 

• the internal masses mi .jv are oriented clockwise as well, with mi between p^ and pi. 
Let K he Si multi- index such that: 

Bk{12) = BK{plml,ml), 
Cif(123) = CK(pl,pl,(pi+P2T,ml,ml,ml), 
Dx(1234) = Dk(pI,pI,pI, (Pi +P2+ P3)^ (pi + P2)^ (P2 + P3)^ mi, ml, ml, ml). 
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The reduction of tensorial functions into scalar functions can then be done according to the 
following standard notations: 

5^(12) = J9^5i(12), 
B'^'^iU) = p^p^B2i(12) + ^'^^522(12), 

= p^'Cii(123)+p^Ci2(123), 

= p5'KC2i(123) + p^p^C22(123) + pS^p^^C23(123) + ^'^'^C24(123), 

= (5'^''p? + y'^X+5>nCooi(123) 

+ (^'^>^ + ^'^>^ + 5>^)Coo2(123) 

+ p^KP?C'ni(123)+p^p^p^C222(123) 

+ {P'M + p'iP2P'i + P^pK)Cii2(123) 

+ (pM +P»^)Ci22(123), 

^''(1234) = p^Dn(1234)+p^Di2(1234)+p^Di3(1234), 

£"^'^(1234) = p^KAi(1234)+p^p^D22(1234)+p^p^D23(1234) 

+ pW/^24(1234) +^{''^3^^^25(1234) +p5'^p^>L'26(1234) +5''^L>27(1234), 

^'^'^^(1234) = ((?^>? + (?>5' + (7>^)Dooi(1234) 

+ (^7''>^ + ^7>^ + ^7>^)I)oo2(1234) 

+ {g'''p's + 9'"'p'3+9'Y3)Doo3{1234) 
+ E Prp;p^A,fc(1234). 

l<ijk<3 

In the reduction of D^'-''', the sum is over all triplets {i,j,k) with repetitions, i.e. 3^ = 27 
terms. By construction, the coefficients Z^jj^ are invariant under index permutations. 

More details about this approach and about the reduction of the PV tensorial integrals into 
scalar ones can be found in [19]. 



C"^(123) 
C'''^(123) 
C^'^''(123) 



B.2 LoopTools definitions 

Sometimes, as for instance in the LoopTools library [16], it is convenient to use another 
notation and to introduce momenta ki given by: 

ki = pi, 

k2 = P1+P2, 

^3 = P1+P2+P3-- 



N 



kN = ^Pi 

1=1 

In that case, the tensorial coefficients are: 

Bi" = k'lB^, 
B^^ = klklB^.+g^'^B^. 
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ij=l,2 

ijl=l,2 i=l,2 



= k^D^ + k^D^ + k^D^, 
D^'' = Y: k^k'^D^j + g^'^D^, 

l<ijk<3 

where C^^, C^^; and D^^ are completely symmetric. 

By expanding these equations and by then comparing all their terms to those arising from 
the reduction of standard tensorial functions, one can find the relations that exist between 
the standard PV functions and those which are computed in the LoopTools library. 



For the 2 point functions, one obtains: 



For the 3 point functions, one obtains: 



Cii 




Cl2 


= c'2 




• 


C21 


— ^11 


C22 


— ^22 


C23 


— ^12 


C2A 


- 

— '^oo 



B21 — -Bn 
B22 — -Boo 



■^2 



+ 2Ci2 + C2 



+ C22 



C'ool — C'mi + C/^ 



001 — <^001 "I" ^^002 

1 

002 — <^002 



Cm — Ciii + 3Cfi2 + 3C^2 + ^222 

C222 — C222 

C112 — C^l2 + 2(7^22 + C222 

C122 — C122 + C222 
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For the 4 point functions, one obtains: 













+ 




• 


^3 


D21 


= 


D^i + Dg + Dg + 2{Dg + Dg + Dg) 


D22 




^22 + 2^23 + ^33 


D23 






D24 




Dg + Dg + Dg + 2D^^ + Dg 


D25 




Dg + Dg + Dg 


D2& 




-'^23 ' ^33 


D27 




-^00 



Dqoi = -DqOI + -^002 + -^003 

D002 = -D002 + -D003 

-D003 = -D003 

Dm = -Dfn + 50^2 + 3-Dfi3 + 30^2 + 6-0^23 + 3-0^33 + D222 + 3-D223 + 3/^233 + -^333 

-D112 = -Dn2 + -0^13 + 2Df22 + 4-0^23 + 2L'f33 + D222 + 3-D223 + 3-D233 + -D333 

-D113 = -Dn3 + 2Dg^ + 2Dg^ + D223 + 2-D233 + -D333 

-D122 = -0^22 + 2-0^23 + -^133 + -^222 + 3-D223 + 3-D233 + -0^33 

-D133 = -0^33 + -D233 + -^333 

-D123 = -0^23 + -^133 + -^223 + 2-D233 + -0^33 
• 

D222 — D222 + 3-D223 + 3-0^33 + -D333 

-D223 = -^223 + 2-D233 + -D333 

-D233 = -D233 + -D333 

-D333 = -D333 
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